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ABSTRACT 



Context. Shells in Elliptical Galaxies are faint, sharp-edged features, believed to provide evidence of a recent (~ 0.5 -2x10' years ago) merger 
event. We analyse the Globular Cluster (GC) systems of six shell elliptical galaxies, to examine the effects of mergers upon the GC formation 
history. 

Aims. We examine the colour distributions, and investigate differences between red and blue globular cluster populations. We present luminos- 
ity functions, spatial distributions and specific frequencies (5,v) at 50 kpc radius for our sample. 

Methods. We present V and I magnitudes for cluster candidates measured with the HST Advanced Camera for Surveys (ACS). Galaxy back- 
ground light is modelled and removed, and magnitudes are measured in 8 pixel (0.4 arcsec) diameter apertures. Background contamination is 
removed using counts from Hubble Deep Field South. 

Results. We find that the colour distributions for NGC 3923 and NGC 5982 have a bimodal form typical of bright ellipticals, with peaks near 
y - / = 0.92 ± 0.04 and V - / = 1 . 1 8 ± 0.06. In NGC 7626, we find in addition a population of abnormally luminous clusters at M; = - 1 2.5 . In 
NGC 2865 we find an unusually blue population, which may also be young. In NGC 1344 and NGC474 the red cluster population is marginally 
detected. The radial surface density profiles are more flattened than the galaxy light in the cores. As already known, in NGC3923, which has a 
high 5ai of 5.6, the radial density distribution is more shallower than the diffuse galaxy light. 

Conclusions. The clusters in NGC 2865 and NGC 7626 provide evidence for formation of a population associated with a recent merger. In the 
other galaxies, the properties of the clusters are similar to those observed in other, non-shell, elliptical galaxies. 



Key words. Galaxies: elliptical and lenticular, cD; Galaxies: star clusters; Galaxies: photometry; Galaxies: interactions 



1. Introduction 

In current galaxy formation models, most ellipticals have 
already formed at z > 2 JEllisetaL 19971 IPeebles 20021 
Ivan D okkum et al. 2004J. It is unclear whether all globular 
cluster (GC) systems were also formed at this early epoch 
or if substantial numbers are still forming today. An impor- 
tant diagnostic is the existence of bimodality in the colour 
distribution of GCs, present in many early type galaxies 
( |Zepf & Ashman 1993| IWhitmore et al. 19951 1. Generally this 
is explained as being due to metallicity differences indicating 
two or more populations of GCs. Several theories have been 
proposed to explain the origin of bimodality: 



* Based on observations made with the NASA/ESA Hubble Space 
Telescope, obtained at the Space Telescope Science Institute, which 
is operated by the Association of Universities for Research in 
Astronomy, Inc., under NASA contract NAS 5-26555. These obser- 
vations are associated with program G09399 



Merger scenarios ( IToomre 19771 ISchweizer 19871 
[Ashman & Zepf 1992^ in which the metal-rich GCs were 
created in gas rich mergers. Since most star formation occurred 
at early epochs, this means in general that the metal-rich GCs 
are also old. However, this scenario also suggests that GCs 
can still be forming today in mergers. This is supported by 
observations of young cluster-like objects in cuiTent mergers 
in action or possible merger remnants like NGC 4038/39 
(IWhit more & Schweizer 19951 IWhitm ore et al. 1999 1, NGC 
3921 (Schw eizer et al. 1996 1, NGC 7252 (Miller et al. 1997>, 
NGC 1316 HGoudfrooij et al. 2001a| jGoudfrooij et al. 2001b| 
IGoudfrooij et a l. 2004> , NGC 1700 and NGC 3610 
(Whitmore et al. 1997 . 



Two others models explaining bimodality are the accretion 
model dCote et al. 19 98 1 and the multiphase formation model 
(Forbes et al. 19971 HaiTis et al. 1998). In both models all GCs 
are old. The first model produces bimodality by accreting and 
mixing metal poor GCs from dwarf galaxies with the more 
metal rich GCs of the massive host galaxy. Cannibalism by our 
own galaxy of the Sagittarius and Canis Major dwarf galax- 
ies and their clusters ilbata et al. 1995IIForbes et al. 2004a> and 
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observations of large numbers of dwarf galaxies around gi- 
ant galaxies are cited as supporting this scenario. The second 
model explains the bimodality as the result of two phases of 
GC formation in the initial collapse and formation of a galaxy. 
The metal poor clusters, and a small proportion of the stars 
form in the initial gravitational collapse, then the metal rich 
clusters and the bulk of the stellar component form from en- 
riched gas in a second collapse phase about one or two Gyr 
later. [Strader et al. (2004)| argue that this "in situ" model of GC 
formation is in better agreement with their observations of the 
correlation of the colours of the metal-poor populations with 
galaxy luminosity. 

A combination of these different scenarios is used in the 
hierarchical merging model of Beasley et al. (2002)| who un- 
dertook semi-analytical simulations of GC formation. In this 
model the metal-poor GCs are old and formed in cold gas 
clumps, the metal -rich ones are created later in merger events. 
In the hierarchical build up of galaxies, accretion of GCs will 
also take place. These simulations are able to reproduce the 
many variations in the colour distributions of GC systems ob- 
served in elliptical galaxies as well as the observed L - N,o, 
relation. 

Recently, Yoon et al. (2006), showed that the apparent bi- 
modality in globular cluster colours not necessarily implies a 
bimodal metallicity distribution. The nonlinear nature of the 
metallicity-to-colour transformation could cause a single old 
population with a unimodal metallicity distribution to look bi- 
modal. This model is attractive because it gives a very simple 
explanation for the observed distributions and could simplify 
theories of elliptical galaxy formation. However, the observa- 
tions of recent GC formation, as mentioned above will likely 
sometimes disturb the predictions made by this model. 

There are many examples of multi-colour and spectro- 
scopic data for the GC systems of 'normal' elliptical galax- 
ies, and these generally give old ages for both blue and 
red populations (M49: IPuzia et al. 19991 ICohen et al. 20031 
NGC 1399: IForbes etal. 2001bl M87 : ICohen et al. 19981 NGC 
1052 and NGC 7332: Fbrbes et al. 2001a and in a sample 
of early-type galaxies: Strader et al. 2005j . This is in contrast 
to studies of the well-known bimodality of the colour dis- 
tribution of the GCs of the LMC f Gascoigne & Kron 1952| 
van den Bergh 1981, 1991) which is primarily an age ef- 
fect, understood in terms to the evolution with time of clus- 
ters in the colour-colour diagram (IFrenk & Fall 1982t . possi- 
bly combined with a relationship between age and metallicity 
(BattinelH & Capuzzo-Dolcetta 1989; IGirardi e t al. 1995 1. In 
ellipticals, the majority of GCs formed at high redshift (z>2.5), 
whichever formation mechanism is dominant. However it is 
important to examine evidence that recent merger events can 
produce enhanced populations of young clusters, as this pro- 
cess may have been much more important in the early universe. 

This study focused on a sample of elliptical galaxies with 
faint stellar shells in their envelopes. Shells and ripples in el- 
liptical galaxies ( |Malin & Carter, 1980..MaUn & Carter, 1983| 
ISchweizer & Seitzer 1992t are faint, sharp edged stellar fea- 
tures in the envelopes of these galaxies which are the remnants 
of the stellar components of minor mergers in the compara- 
tively recent history of the galaxy ( |Quinn 1984^ . Typical dy- 



namical ages of the shell systems are ~ 0.5 - 2 x 10^ years 
JNulsen 19891 [Hemquist & Quinn 1987| . If we can detect a 
population of clusters then the age of this population will pro- 
vide an independent estimate of the age since the merger, as- 
suming that the same event produced both the shells and the 
young clusters. In NGC 2865, one of our sample, Hau et al. 
(1999) find that the age of a nuclear starburst model for the 
young stellar population in the core of this galaxy is much older 
than the dynamical age of the shells, although they do find a 
better coiTespondence between ages for a model involving trun- 
cation of ongoing star formation. In NGC 1316, Goudfrooij 
et al. (2001b) find a cluster system with an age of 3.0 ± 0.5 
Gyr, consistent with the age of the nuclear stellar population. 
The shell systems of both NGC 2865 and NGC 1316 are com- 
plex. Type II or III systems (P rieur 199 0) whereas the dynam- 
ical age estimates are more directly applicable to simple phase 
wrapped. Type I shell systems. 

Two techniques to investigate possible age differ- 
ences between different GC populations involve measur- 
ing radial density profiles, and globular cluster luminos- 
ity functions (GCLFs). Radial density profiles generally 
show a flattening of the globular cluster density pro- 
file near the centre, when compared with the profile 
of the background galaxy light fLa uer & Kormendy 1987} 
|Capuzzo^ olcetta & Donnarumma 2001). Mechanisms which 
could cause a depletion of the cluster population near the cen- 
tre are dynamical friction, which causes clusters to spiral in to- 
wards the centre ( Tremaine et al., 1975 Pesce et al. 1992> . and 
destruction by tidal shocks as the clusters pass close to the nu- 
cleus dOstriker et al. 1989l|Capuzzo-Dolcetta & Tesseri 1997) . 
This second process operates preferentially in triaxial poten- 
tials, and on clusters on radial orbits. If we can identify a popu- 
lation of younger clusters, created during a recent merger, then 
the density profile might extend further into the centre as the 
clusters have had less time to disrupt. However, this does de- 
pend upon the orbital structures to be the same, if one or other 
population were on predominantly radial orbits, then this would 
cause a stronger flattening of the core. 

The same process can lead to evolution of the mass func- 
tion of the clusters and hence the GCLF ( Fall & Re es 19771 
IGnedin & Osrik er 1997 Fafl & Zhang 2001 1. Lower mass 
clusters are preferentially destroyed, leading to the well known 
turnover in the GCLF. A younger population would contain 
more low-mass clusters, and the GCLF should be closer to 
the original mass function, which might be a power law. 
Observations of young cluster systems such as NGC 1316 
( jGoudfrooij et al. 20 01b i and NGC 7252 JMiUer et al. 1997) 
do show power law GCLFs. 

In this paper we analyse the properties of the GC sys- 
tems of six shell elliptical galaxies using optical V and I 
data from the Advanced Camera for Surveys (ACS) on HST. 
Although the observations were optimised to study the shell 
structures rather than the GCs, the magnitudes, colours and 
spatial distribution will indicate whether GC systems of any 
of the shell galaxies differ from those of normal early type 
galaxies. We will also provide data on three GC systems 
which have not been studied before: NGC 474, NGC 1344 and 
NGC 2865. The last galaxy is known as a recent merger rem- 
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nant ( IHau et al. 1999IISchirmnovich 1995> . The GC systems of 
three of our galaxies have been studied before (NGC 3923, 
ground-based ( |Zepfetal. 19941 [Zepf et al. 1995^; NGC 5982 
and NGC 7626, WFPC2 on the HST (Forbes et al. 1996 1. The 
larger field of view and higher resolution of ACS will provide 
more detections and more complete knowledge of these GC 
systems. 

In Section 2 and 3 the observations and data reduction are 
described, which include the detection and selection of the 
GCs, calculation of completeness levels and photometry of the 
globular clusters. Section 4 presents the V-I colour distributions 
and spatial distributions of the globular cluster systems. Section 
5 describes the globular cluster luminosity function which is 
used in Section 6, where the specific frequencies are calculated. 
Finally in Sections 7 and 8 we discuss the results and present 
the conclusions. 

2. Observations and Data Reduction 

The six shell galaxies (see Table 1) were observed with the 
ACS_WFC camera between July 2002 and January 2003 with 
the filters F606W (V) and F814W (I), with CR.SPLIT=2. The 
ACS camera contains two CCDs of 2048 x 4096 pixels, each 
pixel having a size of 0.049" pixel ' resulting in a field of view 
of 202" X 202". Exposure times were on average 1000s, see 
Table 2. The inner 8 pixels of NGC 474 and the inner 24 pixels 
of NGC 2865 were saturated in both V and I. 

Standard reduction was carried out in the IRAF+STSDAS' 
environment, using the packages CALACS and PyDrizzle. 
These are provided by the Space Telescope Science Institute 
(STScI). CALACS processing includes bias and dark subtrac- 
tion, removal of the overscan regions, flat fielding and cosmic 
ray rejection. The default pipeline is not efficient at removing 
cosmic rays when the image is filled by a large galaxy. To solve 
this problem we changed the value of CALACS pipeline pa- 
rameter SCALENSE from 0.3 to 0.0. Setting SCALENSE to 
0.0 increases the probability of removing good stellar data^. 
However, this is only true for empty fields, which is not the 
case for our data. The resulting images were further processed 
by PyDrizzle, which removes the geometric distortion of the 
ACS optical configuration. Finally, after drizzling the images, 
the IRAF package LA.COSMIC ( Ivan Dokkum 2001) was used 
to remove any remaining cosmic rays, which still affected sev- 
eral hundreds of pixels in each image. The standard ACS pho- 
tometric calibration was used ( S irianni et al. 2005j to obtain 
Johnson V and Cousins I magnitudes. The following transfor- 
mation formulae were applied: 

Vj = m(f606W)-H26.331-H0.340*(y-/)yc-0.038*(y-/)yc(l) 

Ic = m(F814W)-H25.496-0.014*(y-/)yc+0.015*(y-/)5c(2) 

Here m(F606W) and m(F814W) are ACS Vega instrumen- 
tal magnitudes and V and I are in Johnson and Cousins systems 

' IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 

^ http://www.stsci.edu/hst/acs/documents/newsletters/stan0301.html 



respectively. The Johnson and Cousins magnitudes were cor- 
rected for galactic extinction using the values of Schlegel et al. 
(1998), which are listed in Table 1. 

3. Data Analysis 

In this Section we describe how the globular cluster source 
catalogues were obtained and describe their characteristics in 
terms of photometric errors and completeness, which are used 
in the further analysis of the data. 

3.1. GALPHOT 

Information about the morphology of the galaxies was ob- 
tained by using the ellipse fitting task GALPHOT (see 
|j0rgensen et al. 1992| l; it returns information such as elliptic- 
ity, position angle, surface brightness and the C3,C4,S3,S4 co- 
efficients C Carter 197 8 1, all as a function of radius. A galaxy 
subtracted residual image is also returned, which we used to 
extract the GCs. 

In the GALPHOT processing, background galaxies, point- 
like objects, dust lanes, bright pronounced shells and additional 
bad data were masked out in an iterative way. Remaining faint 
shell structures, having a brightness typically not more than 5% 
of the galaxy light, do not severely affect the final results. The 
best fits were obtained by allowing the center, position angle 
and ellipticity to be free parameters. In two cases: NGC 2865 
and NGC 5982 the central regions could not be subtracted in a 
proper way, these regions correspond to circles with diameters 
4" for NGC 2865 and 24" for NGC 5982 respectively. Note 
that the inner 0.3" is saturated in the center of the latter galaxy. 

Light profiles were obtained by plotting the surface bright- 
ness for each fitted ellipse as function of radius. The outer 
parts of the light profiles are severely affected by uncertainty 
in the determination of the background. It is diflicult to deter- 
mine a reliable background values from the ACS images them- 
selves, since the galaxies fill the whole field. Fortunately, for 
five galaxies we found optical wide field data in the R band in 
the ESO archive of WFI@2.2m. The WFI camera has a field of 
view of 34' x 33', much larger than the galaxies. 

We used the ASTRO-WISE system ^ 
( |Valen tijn & K uijken 2004| l to reduce the WFI images. 
After subtracting a constant background value from the WFI 
images, GALPHOT was applied to them. The ACS background 
values were determined by matching ACS light profiles to the 
WFI fight profiles. For the galaxy without WFI data, NGC 
5982, we extrapolated from the outer points of the image using 
a |de Vaucouleurs (19 48) r'^"^ law. The calculated values for 
the backgrounds are listed in Table 2. The final GALPHOT 
results describing the morphology of the six galaxies will be 
presented in a future paper (Sikkema et al. in preparation). 

3.2. Globular cluster candidates 

The galaxy subtracted, residual images were used for source 
detection and photometry of globular cluster candidates 

3 www.astro-wise.org/portal 
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Table 1. Properties of six shell galaxies. Data in columns 2-9 from Roberts (1991). 1st column: Galaxy name; 2nd and 3rd 
column: Right Ascension and Declination; 4th and 5th column: galactic longitude and latitude; 6th and 7th column: extinction 
coefficients in V and I magnitudes from Schlegel et al. (1998); 8th column: Morphological type (throughout this paper we 
assumed that NGC 474 is elliptical and not SO CHau et al. 1996'); 9th column: total apparent V magnitude; 10th column: SBF 
distance modulus in I, corrected for extinction. (Tonry et al. 2001 1, f: distance moduli taken from Roberts et al. 1991 determined 
from HI velocity data, corrected for galactic rotation and restframe of the Local Group using i/o=75km/s/Mpc. 11th column: 
distances adopted in this paper using column 10 except $: GCLF distances for NGC 1344 and NGC 3923, calculated in Section 
5.2.; 12th column: Dn-o" distances from Faber et al. (1989); 13th column Central velocity dispersion cr; typical eiTors are 10 
km/s (HYPERLEDA'*) o: d(Dn-cr) for NGC 2865 is probably wrong, since NGC 2865 exhibits a central depression in cr due to 
a rotating disk JHau et al. 1999> . Extrapolating his data gives cr=230km/s (d(Dn-cr)=34.7Mpc), much larger than the value used 
by Faber et al. (cr=168km/s). 
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Table 2. Observational characteristics. 1st column galaxy name, 2st and 3nd column: V and I exposure times in seconds, 4th and 
5th column: 80% completeness levels in V and I, 6th and 7th column: adopted background value in counts 



(GCCs). SExtractor JBertin & Arnouts 199 6^ was used for this 
purpose. The keyword BACK_FILTERSIZE was set to 5; the 
keyword BACK_SIZE was set to the relatively small value of 
24, to account for background variations caused by broad shell 
regions, diffracted light from bright stars, and dusty regions. 

3.3. Photometry 

Because we are analysing point like objects, the magnitudes are 
measured in fixed apertures, and we tested apertures with diam- 
eters of 4,6,8, 10, 12, 14, 16, 1 8,20 pixels. The quality of the pho- 
tometry was analysed by adding objects to the residual images 
using the IRAF task MKOBJECT (the object was extracted 
from an ACS image and is a typical PSF object with a few thou- 
sands counts at its peak). We did this 5 times, each time intro- 
ducing 350 objects with a uniform magnitude distribution and 
positions (generated by IRAF/STARLIST. Catalogues were ex- 
tracted by using SExtractor again with the same parameters 
as before. The resulting catalogues in V and I were then as- 
sociated, keeping only associations within positional error el- 
lipses of lxFWHM=2.4 pixels. Comparing input STARLIST 

http://www-obs.univ-lyonl .fr/hypercat 



data with output data gives information about photometric er- 
rors and completeness levels. 

3.4. Completeness 

We define completeness as the ratio of recovered objects di- 
vided by the original number of added objects measured in a 
particular magnitude bin. FigureQlshows a typical result of this 
procedure for the I band residual image of NGC 3923. On the 
horizontal axis are the input STARLIST magnitudes. The com- 
pleteness is plotted on the vertical axis. The horizontal lines in- 
dicate 50% and 80% completeness levels. The results for each 
image are shown in Table 2. 

3.5. Photometric errors and aperture selection 

The photometric errors were obtained by subtracting in- 
put STARLIST magnitudes from output aperture magnitudes, 
where extreme outliers (mainly due to false associations) in 
each magnitude bin were removed using 5x sigma clipping. 
Typical results, using the NGC 3923 I band residual image and 
an 8 pixel aperture, are shown in Figure |2] The upper panel 
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Fig. 1. Histogram completeness for NGC 3923 I band image, 
using the gc-type object. Vertically: completeness ratio; hori- 
zontal scale: instrumental ACS magnitudes. Dashed lines indi- 
cate 80% and 50% completeness levels. 
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Fig. 2. Photometric errors for artificial pointlike objects in the 
GALPHOT NGC 3923 1 band residual image. Top: RMS errors 
in magnitudes in each bin; bottom: Off-sets from STARLIST 
artificial objects using an aperture of 8 pixels. 



shows RMS eiTors in magnitudes in each bin; off-sets from 
STARLIST magnitudes are plotted in the lower panel. The 
magnitude eiTors depend on the different apertures in various 
ways: 

- The rms-errors at the 80% completeness levels increase for 
larger apertures: from 0.08 mag using a 4 pixel aperture, to 
0.18 mag using a 20 pixel aperture. 

- The measured magnitudes show systematic off-sets for 
each aperture; larger apertures give smaller offsets. We as- 
sume that the offsets of the real data are the same in all of 
our twelve images. 

- The off-sets vary within each magnitude bin; however the 
variations are larger for large apertures. 

Considering these three eiTor sources leads to a optimum 
aperture diameter of 8 pixels, which was used in all remaining 
analysis. The photometric eiTor is k 0.10 mag. near the 80% 
completeness levels and the aperture correction is 0.26 mag. 
Figure 121 shows an example using this aperture for the NGC 
3923 I band residual image. 

Before doing the photometry, we checked whether glob- 
ular clusters are resolved in our galaxies. Table 1 shows that 
our closest galaxies, NGC 1344 and NGC 3923, have a dis- 
tance of about 20 Mpc. At this distance, typical galactic glob- 
ular clusters, which have a half total light diameter of 6 pc 
( |van den Bergh 1994| , will have an angular size of somewhat 
more than 1 ACS pixel. We conclude that most globular clus- 
ters will be unresolved. 



3.6. Selection of GCCs 

We made use of various SExtractor keywords to select the 
globular cluster candidates. We kept all objects which have 
ELONGATION < 1.4 (keep round looking objects), 1.8 < 
FWHMJMAGE < 5 (removing most spurious objects and ex- 
tended objects) and FLAGS > (removing objects which are 
blended, saturated or placed near other bad pixels). Objects 
fainter than the 80% completeness levels were also removed. 
The remaining lists of objects in V and I were associated using 
positional error ellipses of lxFWHM=2.4 pixels, which again 
removed many objects. The final source list was further cleaned 
by visual inspection of the residual images, thereby removing 
any false data, for instance bright galaxy cores, sources located 
near borders or detections in spikes of bright stars. Figure |3] 
shows the distribution of the remaining objects. 

Before using these data, we note that the source catalogue 
is still contaminated by two sources: foreground stars and back- 
ground galaxies. Below, an estimate of these numbers is made. 

1. Galactic foreground stars. The number of foreground stars 
depends upon galactic coordinates. Table 1 shows that 
NGC 2865, with the lowest galactic latitude of 19°, is likely 
to be affected most by foreground stars. We estimated the 
number of stars in each field by using The Besanfon model 
(>Robi n et al. 2003| l. We used the following input for the 
model: 

- An error polynomial - a polynomial was fitted to our 
rms-error curve (top panel in Fig.|2j and given as input. 

- Magnitude limits = our 80% completeness levels 

- Field of view = ACS field size 
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pleteness limits. Figure^shows that we reach 100% complete- 
ness at these brighter magnitude limits. Another advantage of 
doing this, is that the photometric errors are smaller at these 
limiting magnitudes (0.07 magnitude instead of 0.10). The fi- 
nal V-I colour distributions are depicted in the panels of Figure 
|4] The upper part of each panel shows the distribution of objects 
in the colour - absolute magnitude plane. Here, the vertical axis 
is absolute I-band magnitude, which was determined using the 
distances listed in column 11 of Table 1 . The lower part shows 
V-I histograms with binsize 0.075 mag. In these panels the up- 
per, thin-line histogram represents the data, the lower, thick- 
line histogram represents the residuals from the fits described 
in Section l4!2l 

The expected number of stars in the Besan§on catalogues is 
very low compared to the number of GCCs and will not influ- 
ence any conclusion drawn from the figures. In the NGC 2865 
data, the galaxy at the lowest galactic latitude, there is a small 
bump near V-I=1.6, which is also expected from the Besangon 
model. The curves and thick lines are Gaussian fits to the his- 
tograms and their residuals respectively (see Section B3t . 



Fig. 3. Location of globular clusters, using ACS pixel coordi- 
nates (pixelsize = 0.05 arcsec; the open square represents the 
center of the galaxy. 

- Galactic coordinates. 
The model returns catalogues with the expected number of 
stars and their V and I magnitudes. Generally we find that 
the number of stars is negligible compared with the num- 
ber of GCs. A more detailed discussion of the eff'ect of this 
contaminant is given in section 4. 1 . 
2. Unresolved background galaxies. An estimate of the num- 
ber of these contaminants was made by using the Hubble 
Deep Field South (IWilhams et al. 19961 . The images are 
publicly available and were observed in the same pass- 
bands as our data. Since the HDFS images are much deeper 
than our images, we dimmed the HDFS in order to reach 
our 80% completeness levels. We did this by dimming the 
HDFS V and I images with 2.8 and 2.6 magnitudes re- 
spectively, by multiplying with iQ-<'^'"''s>/2.5 ^ restore 
the noise levels of the original images, we added Gaussian 
noise. Next, SExtractor was applied with the same param- 
eters as before, selecting objects using the same selection 
criteria as before and then associating the objects. We find 
15 sources within 0.55 < V - I < 1.45, which are evenly 
spread in V-I. This number, together with the number of 
foreground stars, will be used in Section 4.4 to estimate a 
contaminating background density valid for each galaxy. 

4. V-I distributions and spatial distributions 

4.1. V-I distributions 

To appear in our histograms an object must be detected in both 
passbands, so to avoid colour bias in the selection, we selected 
only objects from our source catalogue which are 0.3 magni- 
tudes brighter in both V and I than the respective 80% com- 



4.2. Components of tlie colour distributions. 

Inspection of the histograms in Figure |4] suggests that the 
colour distributions of GCs for NGC 3923 and NGC 5982 have 
the bimodal form which is normal for bright ellipticals, with 
blue and red peaks near (V - /) = 0.92 and (V - /) = 1.18 re- 
spectively. This bimodal form is parameterised conventionally 
as the sum of two Gaussians. In section we investigate the 
ways in which the other four histograms differ from this form. 

A check for colour bimodality was made by apply- 
ing the KMM algorithm ( A shman etal. 1994> and DIP test 
( ,Gebhardt& Kissler-Patig 1999) 1 on data points between 0.55 < 
V-I < 1.45. KMM was used in standard mode (fitting two 
Gaussians with equal cr's). It returns the value P, which indi- 
cates if a distribution is better characterised by a sum of the 
two Gaussians than a single Gaussian. Table 3 lists the output 
which consists of peak values, number counts in each Gaussian 
and the P-value. 

The KMM analysis shows that all of the colour distribu- 
tions are better fit by a double than a single Gaussian, but in 
NGC 1344 and NGC 474 the numbers of clusters in the red 
(metal-rich) Gaussian are too small to be statistically signifi- 
cant, so we plot in Figure |3 only the blue Gaussian for these 
two galaxies. For NGC 7626 and NGC 2865 neither a single 
nor double Gaussian provides a good fit to the histogram, and 
the structure is more complex. However we do show in Figure 
0]the blue and red Gaussians generated by the KMM algorithm. 

The obvious colour bimodality for NGC 3923 
and NGC 5982 is confirmed by using the DIP test 
(|Gebhardt & Kissler-Patig 1999| l. This test calculates the 
probability of a dip occurring in a supposedly bimodal distri- 
bution. Applying this test to the colour distributions of our six 
galaxies gives significant bimodality for NGC 3923 and NGC 
5982 with the dip probabiHty Pji,, = 0.99 and Pji,, = 0.92 
respectively. 
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Fig. 4. Histograms and colour diagram of V-I colours of the globular cluster candidates. All histograms are better fit by two 
Gaussians than one, but for NGC 1344 and NGC 474 we plot only the blue component as a Gaussian fit, as the numbers 
in the red peak are not significant. For the other four galaxies we plot the two Gaussians fit by the KMM algorithm. Residual 
histograms (thick lines) represent the difference between the data and the sum of the Gaussian fits. Vertical, dashed lines represent 
the separation of the data into red and blue groups at (V - /) = 1 .05. 



NGC 3923 and NGC 5982 have been studied before. 
NGC 3923 has been observed from the ground in the 
Washington system ( |Zepf et al. 1995^ . Like us, they found a 
bimodal distribution, which, converting their data into V-I 
lIForbes & Forte 200 It and correcting to our extinction scale, 
give colours of 1.01 (+0.05) and 1.21 (+0.05), an offset 
of +0.06 with respect to our data. Data on NGC 5982 
IIKundu & Whitmore 200 lai also revealed a bimodal distri- 
bution; applying our extinction scale to their data gives V-I 
colours of 0.96 (±0.03) and 1.15 (±0.03). While the blue peak 



is the same as ours, their red peak is 0. 1 magnitude bluer. We 
attribute this to our much larger GC sample. 

In the further analysis we distinguish between red and blue 
sample by cutting the distribution at (V-I) - 1.05, in order 
to analyse whether these two populations have different prop- 
erties. 

If we fit the V-I histogram of NGC 2865 with 2 Gaussians, 
as we do for NGC 3923, NGC 5982 and NGC 7626, the blue 
peak extends much further to the blue than in other galaxies. 
We suggest in section0that this can be attributed due to a pop- 
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NGC 474 
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1.22 
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0.035 


1.2-2.2 










-1.50 


0.14 


-1.87 


0.03 


NGC 1344 


0.92 


1.21 


268 


49 


0.004 


1.5-2.2 










-2.02 


0.17 


-2.23 


0.02 


NGC 2865 


0.85 


1.12 


192 


111 


0.039 


1.35-2.15 


-3.05 


0.24 


-2.04 


0.17 


-2.28 


0.14 


-2.04 


0.01 


NGC 3923 


0.94 


1.16 


253 


390 


0.036 


1.5-2.1 


-0.87 


0.15 


-0.86 


0.15 


-0.90 


0.10 


-1.59 


0.03 


NGC 5982 


0.96 


1.24 


240 


221 


0.000 


1.4-2.1 


-1.98 


0.18 


-1.22 


0.19 


-1.68 


0.13 


-2.23 


0.04 


NGC 7626 


0.92 


1.13 


455 


448 


0.039 


1.5-2.2 


-1.99 


0.14 


-1.42 


0.13 


-1.67 


0.09 


-2.16 


0.04 



Table 3. Left Output of KMM algorithm (explanation: see Section 4.2) with blue and red V-I peaks (columns 2, 3), KMM 
assignments of GCs to blue and red groups (columns 4, 5) and P-value (column 6), See also Fig. 4. Right: Slopes a and errors of 
the surface densities of red, blue (based on a split at V-I=1.05) and all GCs (columns 8-13), determined applying a weighted least 
squares method on ranges listed in column 7. Last column: slope of galaxy light. Surface densities of all, red and blue groups are 
shown in Figures 8, 9 and 10. 




Fig. 5. Globular clusters in V-I (vertical scale) as a function of 
radius. 



ulation of young blue clusters, overlapping in colour with the 
normal old metal-poor population. 

For NGC 7626 KMM gives maxima at 0.92 and 1.13; the 
bimodal nature was already observed by earlier, but less deep, 
observations ( Kundu & Whitmo re 200 la> . However inspection 
and the output of KMM show that the simple bimodal form is 
not a good fit, instead the histogram shows a broad flat peak. 
The middle of this peak ( (V-I) ~ 1 .0) is filled by exceptionally 
bright (Ml < -11) and probably young GCs (see Section FOl . 
The colour magnitude diagram shows that these are distributed 
in several small clumps in, which cannot be explained by ran- 
dom effects: the errors in V-I at these magnitudes (my ~ 23.0) 
are at most 0.04 (Figure |2j and smaller than the distances be- 
tween the clumps. This galaxy shows that the distribution of 
globular cluster colours is not necessary well represented by 2 
Gaussians. 



Fig. 6. Average colours of globular clusters (circles: all GCs) 
and galaxy colour (curve through triangles) as a function of 
radius. NGC 2865, NGC 3923, NGC 5982 and NGC 7626 also 
show data points for blue and red groups. 

4.3. Spatial distributions of the globular clusters 

Figure|3shows the spatial distribution of the GCCs. The cen- 
ters of the galaxies are represented by an open square. The dis- 
tribution of GCs roughly follows the ellipticities of the under- 
lying galaxies. The lack of GCs in the central regions is due to 
dust, bad fits by GALPHOT (see Section|2j or increased noise. 
Figure Is] shows V-I colours of individual GCCs as a function 
of radius. These colours as a function of radius are plotted in 
bins in Figure |6l Here, triangles, represent the galaxy colour 
returned by GALPHOT, and circles the average colours of all 
GCs. For those galaxies for which we can separate the popula- 
tion into red and blue groups (Section we plot as crosses 
the mean colours of the red and blue groups respectively. 

Except for NGC 1344, the average colours of all GCs 
tend to be bluer at large radii, confirming earlier results 
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Fig. 7. Ratio of red vs. blue clusters as a function of radius; 
NGC 5982 and NGC 7626 have strong gradients. 
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Fig. 8. Surface density of all globular clusters, starlike symbols 
show galaxy surface brightness (arbitrary scale). The slope of 
the GC surface density was calculated between the dashed ver- 
tical Unes. 



JForbesetal. 1996\ and reflecting the colour gradients in the 
stellar halos. Where there is a significant red subgroup the clus- 
ter colours in that group match quite well the galaxy colours, 
except in NGC 3923 where they are bluer, as is the case with 
NGC 1052 (Forbes et al. 2001a> . In NGC 5982 and NGC 7626 
the strong blue-ward gradient in the mean cluster colour is 



caused by a gradient in the relative fractions of the red and 
blue groups, as illustrated by Figure where we plot the ratio 
of the red to the blue population against radius. 

4.4. Globular cluster surface densities 

The GC surface densities were calculated in elliptical annuli; 
their ellipticity was the average value in the outer regions as 
determined by GALPHOT. Dusty regions, bad pixels and other 
bad regions were not used when calculating the effective area 
of each annulus, from which the surface densities are derived. 
We accounted for contaminating sources by subtracting a back- 
ground density, calculated as described in Section l3T6l These 
background densities were between 4.1 to 8.6 x 10""^ arcsec"^ 
for the six galaxies. In Figure|S]we plot the radial surface den- 
sity of GCs as a function of radius, the galaxy surface bright- 
ness distribution is also shown in this plot. In Figures |5] and 
[TUl we plot the radial density distribution for the clusters in the 
red group and blue group respectively, for those galaxies with 
a significant red population. 

For all galaxies we note a deficit of clusters in the inner 
regions, with respect to the background light surface bright- 
ness distribution. Possible physical causes of this are discussed 
in Section l73l but here we consider the possibility of greater 
incompleteness in the inner regions causing this deficit. Two 
effects could contribute to greater incompleteness in the inner 
regions: confusion due to crowding; and the increased photon 
noise level in the higher surface brightness regions. 

To test the effect of confusion a low density region of 1000 
X 1000 pixels, containing 43 clusters, was cut from the residual 
image of NGC 7626. As in Section l331 we introduced artificial 
objects at various number densities, using uniform and power- 
law surface density distributions. After using exactly the same 
detection methods and selection criteria as described in Section 
13.61 we found that we start losing objects due to confusion ef- 
fects at densities of 0.15 arcsec"^, where 5% are missed. From 
Figure |8] we see that even the highest density of our galaxy 
sample (i.e. NGC 7626) is still below this value. 

The effect of the galaxy surface brightness on complete- 
ness is due to the increased photon noise in high surface bright- 
ness regions. Within the central 10 arcseconds (13.5 arcseconds 
for NGC 5982) this leads to increased incompleteness and we 
do not plot points within these regions. Outside this, the back- 
ground photon noise from the galaxy is at most 2x the back- 
ground noise, except for NGC 3923 were the photon noise of 
the inner data point is about 4x the background photon noise. 
After again introducing artificial objects and using the same se- 
lection criteria as before, we find that this effect is significant 
only for NGC 3923: the inner two points in Figure |8] should 
probably lie somewhat higher. 

The radial density distributions of most of the GC sys- 
tems follow the surface brightness distributions of the galaxy 
light in the outer parts, but show a deficit at small radii which 
we argue is not due to incompleteness due to confusion or 
photon noise from the galaxy. This is typical of ellipticals in 
general (e.g. [Lauer & Kormendy 1987) IGrillmair et al. 19861 
IPuzia et al. 20041 IForbes et al. 200 lal ISchweizer et al. 19961 
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Fig. 9. Surface density of red globular clusters; starlike sym- 
bols show the galaxy light (arbitrary scale). The slope of the 
GC surface density was calculated between the dashed vertical 
lines. 

and IBrown et al. 2 000 1. The deficit at the centre is often in- 
terpreted as evidence for tidal disruption of clusters passing 
through the core of the galaxy JFall & Rees 1977j . but to be ef- 
fective this process requires the clusters to be on predominantly 
radial orbits JGrillmair et al. 1986llOstriker et al. 19891 

NGC 3923 shows a different behaviour, the surface den- 
sity profile is much shallower than the light profile at all radii 
(already noted by 'Zep f et al. 1994j l. This is typical of the GC 
systems of some brightest cluster galaxies, e.g. NGC 4874 
JHarris et al. 2000 1. However it is unusual for a galaxy such 
as NGC 3923, which is the brightest in a small group. 

Comparison of Figures|9]and^|shows that there are some- 
times differences between the surface density profiles of red 
and blue GCs. To quantify these differences we fit the linear 
relation logicr) - b + ax rhy applying a linear weighted least 
squares fit to the outer points of the density curve. These points 
were chosen to lie between an inner and outer radius. The inner 
radius was defined as the point where the flattening stops, esti- 
mated by eye. The outer radius excludes unreliable data points 
at larger radii with very low number statistics (usually contain- 
ing only 1 or 2 objects in the partial elliptical ring located in 
the very outer corners of the image). We indicate the inner and 
outer radius in the Figures as vertical dotted lines. 

Using the same method and radii we also fitted the galaxy 
light profile. In columns 7 to 15 of table 3 we list the radii 
used and the results of the fits, including errors on the slopes 
(we used the standard recipe and error formulae, i.e. equations 
15.2.6 and 15.2.9 respectively as listed in lPress et al. 19921 

In NGC 7626, the slope of the red GCs is significantly 
steeper than the slope of the blue GCs. Similar but less sig- 
nificant differences are visible in NGC 2865 and NGC 5982. 



Fig. 10. Surface density of blue globular clusters; starlike sym- 
bols show the galaxy light (arbitrary scale). The slope of the 
GC surface density was calculated between the dashed vertical 
hnes. 

We checked if these differences are due to the choice of our in- 
ner radius by shifting the inner radius one data point to the left 
as well as to the right; we found no significant change in slope 
differences. The slopes calculated for the set of all GCs will be 
used in Section|6j where the specific frequency is calculated. 

5. Globular cluster luminosity function in I 

In this section, the GCLFs in the I band are determined. 
If the observations are deep enough to cover the abso- 
lute turnover magnitude (TOM) of M/ - -8.46 + 0.03 
dKundu & Whitmore 2001 at . the GCLF can be used as a dis- 
tance estimator. The GCLF can also be used to estimate the 
number of globular clusters in a system; this number then 
determines the specific frequency S n (Section |6). In those 
galaxies with a significant red population, we also compare 
the GCLF of the red and blue samples, which could be dif- 
ferent, for instance if the samples have large M/L differences 
( Whitmor e et al. 2002> . 

5.1. Determination 

We calculated the GCLF in the I band, which is less affected by 
extinction. The GCLFs were constructed from histograms with 
a binsize of 0.25 magnitudes, approximately twice the photo- 
metric eiTor of the faintest objects. The Besangon catalogues 
were used to correct for contamination by foreground stars. 
The foreground stars (see Section lT6i were cumulatively sub- 
tracted: i.e. if a bin contains only three sources, while there are 
an expected number of five stars, we subtracted the remaining 
two stars from the next bin. We only used objects within the 
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Fig. 11. Globular cluster luminosity functions in the 1 band. 
Dashed and dotted histograms represent the GCLF of the red 
and blue groups respectively. Gaussian fits are also drawn. Only 
the GCLFs of NGC 3923 and NGC 1344 extend past the TOM. 



80% completeness levels and also corrected for incomplete- 
ness. The blue and red GCLFs were calculated used the red 
and blue groups defined in Section l4!2l The GCLFs are shown 
in Figure ^2 where they are divided in blue and red groups, 
we find no significant differences between the GCLFs of red 
and blue populations. The peak in the GCLF of NGC 1344 at 
1=24.5 is due to red objects, caused either by photometric er- 
rors near the faint limit, or more likely a local excess in the 
background. 

5.2. GCLF as a distance estimator for NGC 1344 arid 
NGC 3923 

If the events which formed the shells also formed large num- 
bers of new GCs, then we might find a difference between the 
TOM of the GCLF, firstly between red and blue populations, 
and secondly when compared with the canonical value for nor- 
mal ellipticals of M/^^^_,, = -8.46 (IKundu & Whitmore 200Tal . 
However we only cover this absolute magnitude for the closest 
two galaxies, NGC 1344 and NGC 3923. Using the IRAF tool 
NGAUSSFIT we find for NGC 3923 cr = 1.4 + 0.1, TOM, = 
23.04 + 0.06, and, ignoring the spike in the GCLF at faint mag- 
nitude, for NGC 1344, cr = 1.35 + 0.10, TOM, = 23.26 + 0.10. 

Comparing our GCLF distance estimates for NGC 1344 
(m-M=31.72) and NGC 3923 (m-M=31.50) with the SBF dis- 
tance moduli (from column 10 in Table 1: 31.48 + 0.30 and 
31.80 + 0.28 respectively), shows that the differences of +0.24 
and -0.30 do not exceed the typical errors between these two 
methods dRichtler 2003i) . A PNLF distance estimate for NGC 
1344 dTeodorescu et al. 20051 with m-M ^ 31.40 + 0.18 is in 



excellent agreement with our value. In the remainder we will 
use our calculated distance moduli for NGC 3923 and NGC 
1344. 

5.3. Total numbers of GCs 

A Gaussian fit to the GCLF of the other four galaxies is 
also necessary to estimate the specific frequency (Section |6}. 
Because the magnitude limit is below the TOM of the GCLF, 
we fit only for cr and amplitude, keeping the distance modu- 
lus fixed within the error bars (0.26 mag.). Doing this for NGC 
5982 and NGC 7626 give cr of 1.5 and 1 .25 respectively, and 
good fits. For NGC 474 and NGC 2865 we find cr ~ 1.05, 
which is unusually small, and conclude that the calculated S ^ 
will be uncertain, and that deeper data are required. 

6. Total number of globular clusters and specific 
frequencies 

In this section we estimate the total number of clusters and spe- 
cific frequency S which is a measure of the number of GCs 
per unit galaxy luminosity. It is a strong function of galaxy 
type and environment: late type galaxies and isolated galax- 
ies usually have lower S n than early type and cluster galaxies 
(Harris 1991J. The calculations are straightforward using the 
results of previous sections. 

The specific frequency is defined by 
Parris & van den BerghT98T) : 



5«=M„,10' 



,0.4(Mv+15) 



(3) 



where A^,,,, is the total number of GCs, and My is the ab- 
solute V magnitude of the galaxy, derived from the apparent 
magnitudes and adopted distances (Table 1). 

The total number of GCs is obtained by using: 



^rot — TZ i]^ missed "I" found) 



found 



(4) 



where N found is the detected number of GCs, corrected for 
incompleteness, Ngauss is the expected number of GCs assum- 
ing a Gaussian distribution: 



N, 



1 



gauss 



binsize 



(5) 



with binsize=0.25; amplitude and cr are from the Gaussian fits 
from the previous section . 

Nmissed is the number of GCs expected to be missed due to 
the fact that the galaxy is much larger than the ACS field of 
view. This number is calculated by extrapolating and integrat- 
ing the GC surface density as determined in Section from 
the outer points (listed in column 7 of Table 3), to the point 
where ^4ab reaches 50 kpc and lOOkpc. Error sources used in 
the calculation of S m are: 

1) A photometric measurement error of 0.15 mag. 

2) A distance error of 15%. This error is taken from 
paber et al. (1989)| and represents the scatter in the Dn — cr dis- 
tance estimator. Comparing our adopted distances with Dn — cr 
distances (see columns 11 and 12 in Table 1) shows that this 
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is a reasonable assumption. The distance error was applied to 
both N,niss and absolute V. 

3) The errors in the slope of the GC surface density, listed 
in column 13 of Table 3. 

A possible large error is the unknown value of the slope 
of the GC surface density outside the ACS field of view. We 
assumed this slope to be equal to the fitted GC density profile 
in the outer part of the galaxies. Especially for our two most 
nearby galaxies this assumption is uncertain; these galaxies 
are only covered by the ACS to 15-17 kpc, while we extrap- 
olate until 50 kpc and even 100 kpc. Data on the slope in the 
outer regions of our most distant galaxy, NGC 7626, and sur- 
face density profiles in the literature, support this assumption 
(e.g. Rhode & Zepf 2004, Zepf et al. 199 4|IHarris etar2 000l. 
Finally, for NGC 3923, wide field data ( |Zepfetal. 1994| l show 
a constant slope to a radius of 5.6 arcmin ( 34 kpc). Applying 
a least squares fit to their inner data, covering the ACS field- 
size, gives a slope of -0.82 ± 0.34, comparable to our value 
(-0.90 ± 0.10). Their outer data gives a slightly steeper slope 
of roughly -1.14 + 0.12, which we used to calculate the missed 
GCs in Equation 4. 

Assuming our assumption is justified, we adopt the final 
values for Sn, which are listed in columns 6, 7 (50kpc) and 
8, 9 (lOOkpc) of Table 4. There sometimes are large variations 
between S n calculated at 50 kpc and 100 kpc, this especially 
true for NGC 3923 and is due to the shallow slope of the GC 
density distribution compared with the luminosity distribution. 
[Harris et al. (1 998) also noticed such variations and found a ra- 
tio S°^""'"/Sf = 1.3, with S"^""'" the total Sn and the Sn 
within 40kpc. It is important to keep these different definitions 
of 5 in mind if one compares Sn with other authors. Some 
authors calculate local Sn, others extrapolate to various values 
ranging between 25kpc and 200kpc. 

Integrated out to 50kpc, we find values of S n typical for 
isolated elliptical galaxies, except for NGC 3923 and NGC 
7626, for which we find Sn = 5.6 and = 3.9 respec- 
tively. These values are more typical of cluster ellipticals, in- 
deed NGC 3923 has the highest Sn of any isolated elliptical 
( |Zepf et al. 1994^ . Both of these galaxies are in groups: NGC 
3923 is a dominant group galaxy embedded in an X-ray en- 
velope (Buot e&Canizares 1 998 Pellegrini 1999 1; and NGC 
7626 is the second brightest member of the Pegasus group and 
is also detected in X-rays cO' Sullivan et al. 2001j . 

7. Discussion 

In this section a short summary of the properties of the shell 
galaxies is given; a comparison with the GC systems of 'nor- 
mal' early type galaxies is made; the data are compared with 
predictions of the hierarchical merger scenario simulations 
(Beasley et al. 2002); and possible signs for recent GC forma- 
tion and ages are discussed. 

7.1. The shell galaxies 

All of our galaxies are located in low density re- 
gions; early type galaxies residing in groups and clus- 
ters have a much lower probability of exhibiting shells. 



dMaUn & Carter, 1983| IColbert et al. 200 U . As well as the 
presence of the shells, all six galaxies also show visible 
dust patches and/or lanes, mostly in the central regions. It 
has been shown that all galaxies, except for NGC 1344, 
contain a KDC or show otherwise peculiar kinematic be- 
haviour (Hau et al. 1996" 'Hau et al. 1999" "Cart er et al. T998I 
Emsellem et al. 2004 Balcells & Carter 199 3 1. NGC 474 , 
NGC 5982 and NGC 7626 may be LINERS JHoetal. 19971 . 
Except for NGC 474, unresolved X-ray data are available for 
all galaxies (O' Sullivan et al. 2001 1 and even a 2D X-ray map 
for NGC 3923 ( Buote & Canizares 1998; Pellegrini 1999). In 
this paper we assume that NGC 474 and 1344 are E-type galax- 
ies, although there is evidence in each case that they night be 
classified as SO galaxies. 

The environmental and other properties of the galaxies 
are described briefly below (LGG group numbers from 
ICiarcia 1993> : 

NGC 474 - brightest galaxy of LGG20 (4 members). 
The galaxy is connected with the small spiral NGC 470 via a 
HI tidal bridge dSchiminovich 1997t . Often classified as SO; 
shows some rotation dEmsellem et al. 2004t . May be a LINER. 
NGC 1344 - located at the outskirts of the Fornax Cluster, 
4.9° from the central cluster galaxy NGC 1399. At this 
distance the density is comparable to the field galaxy density 
(Kambas et al. 2000l. Sometimes classified as SO; shows 
rotation within 2 effective radii (ITeodorescu et al. 2005t 
NGC 2865 - isolated galaxy fReda et al. 20041 . remams 
of a rotating HI disk ( Schiminovich 1995}. IHau et a l. 1 9991 
found a KDC and evidence for a young (0.4-1.7 Gyr) stellar 
population; two possible explanations were given: a starburst 
or a truncation of the star formation. 
NGC 3923 - brightest galaxy of LGG255 (5 members) 
NGC 5982 - brightest galaxy of LGG402 (4 members), may 
be a LINER. 

NGC 7626 - second brightest member of the Pegasus group, 
LGG473, of at least 15 members. This probable LINER 
has a radio-jet, directed NE, and a small HI cloud between 
1.5' and 3.0' WSW of the center; the galaxy does not 
contain any HI tidal features JHibbard & Sansom 2003> . 
The core shows orthogonal kinematics to the main body 
(Balcells & Carter 19931 1. with no emission lines, nor signs of 
nuclear young populations. A dust lane in direction ENE is 
visible in the inner 15 ACS pixels of the core. 



7.2. Comparison of the GC systems with normal 
ellipticals 

Our sample is selected on morphological grounds to have un- 
dergone a recent minor merger In this section we investigate 
the effect of the merger event upon the GC systems. All six 
galaxies have a peak at 0.94 + 0.04 in their V-I histograms, in 
common with 'normal' ellipticals. In four of our six galaxies 
we see a distribution of galaxies red-wards of this, this second 
peak again occurs in the majority of normal ellipticals. 
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Galaxy 


^ found 


^ missed 


error 


N 


error 


A',,,, 


error 


S NSO 


error 


S A-lOO 


error 


Ri„ - Rom (kpc) 


My (mag) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


NGC 474 


284 


143 


29 


405 


56 


609 


94 


2.1 


0.5 


2.7 


1.0 


2.5-25.0 


-21.2 


NGC 1344 


352 


93 


25 


367 


28 


464 


44 


1.4 


0.3 


1.5 


0.4 


3.4-17.0 


-21.3 


NGC 2865 


350 


35 


6 


647 


98 


712 


108 


1.6 


0.4 


1.7 


0.4 


4.1-25.9 


-21.6 


NGC 3923 


660 


1254 


162 


860 


67 


2494 


286 


5.6 


1.3 


8.3 


3.6 


3.1-15.4 


-21.6 


NGC 5982 


505 


116 


18 


1218 


161 


1497 


202 


2.6 


0.6 


3.0 


0.9 


5.1-25.6 


-21.9 


NGC 7626 


1051 


141 


16 


2498 


263 


2833 


300 


3.9 


0.9 


4.8 


1.4 


7.4-37.0 


-22.2 



Table 4. Column 2: Detected GCs, Columns 3,4: missed GCs + error; Columns 5,6: number of GCs + error using the values for 
the Gaussians determined in section ls!2l Columns 7,8: total GCs + error applying equation (4); Columns 9,10: Sn within 50 kpc 
with error; Columns 1 1,12: 5 within 100 kpc with error; Columns 13: ACS range where slope has been determined 



All galaxies show a flattening in the GC density profile 
near the central regions. This feature is also a generally seen 
in other GC systems and is attributed to disruption processes 
dOstriker etal. 1989llPesce et al. 1992> . 

The properties of flie GC systems of NGC 3923 and NGC 
5982 are very similar. Both galaxies have blue and red popula- 
tions of comparable size, which peak at roughly the same V-1 
values near 0.95 and 1.20. Using the colours of the GCCs and 
assuming that the red peak has solar metallicity, the evolution 
models of pritze-v. Alvensleben (2()04 ) (to get the metallicities 
of the blue clusters) and Fig. 12 from Goudfrooij et al. 2001^ 
indicate that these systems are old, > 5 Gyr, and evolved sys- 
tems. 

NGC 3923 has been extensively studied by 
|Zepf et al. 1995| They already noted the very high Sn 
for a galaxy located in a low density environment, which we 
confirm. McLaughlin (1999) proposes that high Sn galax- 
ies are probably best explained by taking into account the 
presence of a extended massive X-ray halo. However, NGC 
3923 is probably in contradiction with these results. While 
NGC 3923 has a very shallow GC surface density profile 
and also contains a X-ray halo JBuote & Canizares 19981 
[Pellegrin i 1999t , its X-ray luminosity is about a factor 10 
lower (|0' Sulliva n et al. 200 1 Fukazawaet al. 20061 than the 
galaxies studied by McLaughUn: M49, M87, and NGC1399. 

The colour distributions of NGC 474 and NGC 1344 
(Figure 0} are similar, and appear either to be unimodal, or to 
have a very low red (metal-rich) component. Although KMM 
returns a combination of a blue and red peak as a good descrip- 
tion of the distribution, there are very few GCs in the red peak. 
Such colour distributions are very rare in luminous ellipticals 
(iKundu & Whitmore 200 1 a i. Peng et al. (2006) show that, in 
early-type Virgo galaxies, the red peak is much more promi- 
nent in the more luminous galaxies. However unimodal, blue, 
histograms only become common for galaxies fainter than My 
= -18 (Figures 4 and 6 of Peng et al.), whereas the absolute 
magnitudes of NGC 474 and NGC 1344 are My = -21.17 and - 
21 .07 respectively. They are however the lowest luminosity and 
lowest velocity dispersion galaxies in our sample. 

Environmental differences in the early history of the galaxy 
might lead to differences in the GC (V-I) histograms. In 
the merger model, fewer early mergers would mean less red 
peak clusters. Fewer early mergers might also cause a galaxy 
to retain more of its angular momentum, and thus to be 



an SO or a rapidly rotating elliptical. NGC 474 and NGC 
1344 both have significant rotation (Emsellem et al. 2004 and 
Teodorescu et al. 2005 respectively), and both are sometimes 
classified as SO galaxies. The most luminous galaxies in the 
Peng et al. (2006) Virgo sample with unimodal blue his- 
tograms are NGC 4660, an elliptical with significant rota- 
tion ( Emsellem et al. 2004t : and NGC 4340, an SBO. However 
there are counterexamples in the Virgo sample, such as the 
rotating ellipticals NGC 4564 and NGC 4697 which have 
red peak dominated colour histograms. Insufficient galaxies 
have been studied in sufficient depth to analyse any possi- 
ble correlation between the rotational properties of ellipti- 
cals and their GC colour histograms. Kundu & Whitmore 
(2001b) find a number of unimodal, blue colour histograms 
among a sample of SO galaxies studied with WFPC2 (e.g. 
NGC 2768), but their sample sizes are small. They do how- 
ever find that a lower proportion of their sample of SOs are sig- 
nificantly bimodal, than their equivalent sample of ellipticals 
(IKundu & Whitmore 200 Tat . 



Finally we compare these objects with some prop- 
erties predicted by the hierarchical merger model of 
[Beasley et al. (2002)| In these semi-analytical simulations, 
metal-poor GCs are old and formed in cold gas clumps, while 
the metal-rich are created in time during merger events. Of 
course, accretion of GCs also takes place during the hierarchi- 
cal build-up of galaxies. These simulations are able to repro- 
duce the many variations in the colour distributions of GC sys- 
tems observed in elliptical galaxies. For instance, in Figure 13 
of Beasley et al., a blue -peaked V-I distribution is shown which 
looks similar to our large single peak V-I distributions of NGC 
474 and NGC 1344. His model also roughly reproduces the 
observed Ly - Ntot relation. We plot this relation in the lower 
part of Figure[21together with our six data points. Our systems 
seem to have somewhat fewer globular clusters than predicted. 
This is partly explained by the fact that Beasley et al. do not 
take into account disruption processes, which might reduce the 
number of clusters by 10%-20%. The upper part of the same 
Figure shows the logarithm of the ratio between the number of 
blue and red clusters. NGC 474 and NGC 1344, the two points 
with lowest luminosity, lie somewhat outside the point cloud 
due to their low number of red clusters. 
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NGC 474 
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NGC 1344 
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NP (5.5) 


0.5 
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NGC 2865 


Old plus 0.1 - 1.7t 
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Old plus 0.5 - 1 


NGC 3923 


2 6+°-^ 






NP (10.3) 


0.8- 1.3 


Old system 


NGC 5982 


^^■-'-2.0 


6.8 ± 


1.5 


YP (6.8) 


0.2 
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NGC 7626 


13.9!^;^ 


8.0 ± 


1.3 


NP(1.4) 




Old plus 2-5 



Table 5. Age estimates. Column 2: SSP ages from De nicolo et al. (2005) except t : lHau et al. 1999l and ± IKuntschner et al. 20021 
Column 3: Ages using fine-structure index S method fSchweizer & Seitzer 1992> . Column 4: Old (NP) and Young (YP) systems 
together with their S2 index according to Michard & Prigniel (2004). Column 5: Dynamical ages of type I shell systems using 
Nulsen (1989). Column 6: Approximate age indications from GC (V-I) distribution (this paper) 
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Fig. 12. Top: Ratio of blue and red clusters vs galaxy lumi- 
nosity. Simulated data ( Beasley et al. (2002)| l and the six shell 
galaxies (NGC numbers). Bottom: Total number of clusters 
vs. luminosity for simulated data and six shell galaxies (large 
squares) 



7.3. Possible evidence for recent GC formation in 
NGC 7626 and NGC 2865 

In this section we examine the possibility that the complex (V- 
I) histograms of NGC 7626 and NGC 2865 provide evidence 
for recent GC formation, in addition to the old GC popula- 
tion which gives both the red and blue peaks. In NGC 7626, 
Figure 01 shows that the brightest GCs of NGC 7626, with 
M[ » -12.5 ± 0.2 and V-I=l.l, are more than 1.0 magnitude 
over-luminous with respect to the brightest GCs of the other 
five galaxies (all at about -11.0), and with respect to its own 
universal blue population. It is the most luminous galaxy in our 
sample(Mv - -22.16), and as such the GCLF would be ex- 
pected to extend to the brightest magnitudes, simply because of 



the larger population at the sparsely populated faint end of the 
GCLF. We compare the GC population of NGC 7626 with that 
of NGC 4472 ( [Rhode & Zepf 2001 > , a more luminous galaxy 
(Mv = -22.7) in the Virgo cluster. NGC 4472 has its bright- 
est clusters in the blue (universal) peak, at R=19 (equivalent to 
Mi - -12.5). The brightest clusters in NGC 7626 are about 
0.2 mag brighter and 0.2 mag redder in (V-I) than this, in a less 
luminous galaxy. This group of bright clusters appears in the 
GCLF of NGC 7626 (FigurefTH as a small excess near 7 = 21. 

We can compare this population with that found by 
Whitmore et al. (1997) in NGC 3610, a dynamically young 
elliptical. Their Figure 15 illustrates the evolution of a young, 
metal-rich population compared with an old, metal-poor popu- 
lation. In time, the young, metal rich, GC population will fade 
in luminosity and will become redder. After three Gyr, this pop- 
ulation has become redder than the old population but still has 
several GCs which are brighter than the brightest old metal 
poor GC. This is exactly what is visible in our colour mag- 
nitude diagram. Whitmore et al. quantify the age differences 
using the A{V - 1) vs. AVio diagram (see their Figure 18). This 
diagram represents an age sequence by plotting vertically the 
difference in V-I between the peaks and horizontally the mag- 
nitude difference between the 10th brightest globular cluster in 
the young and old populations. In NGC 7626 AVw is extremely 
difficult to estimate, because the young population is domi- 
nated in number not just by the old metal-poor population, but 
by an old, red, metal-rich population as well. Replacing AVio 
by an estimate of the difference between the magnitude of the 
brightest clusters and the brightest in the blue peak, we esti- 
mate a very tentative merger age of 2 - 5 Gyr. We conclude that 
NGC 7626 appears to have some young GCs, probably created 
in a recent (2-5 Gyr old) minor merger. These young GCs are 
superimposed upon a much richer, bimodal, old cluster popu- 
lation. 

For all galaxies with a substantial red population, we find 
significantly steeper slopes for the red clusters. This is re- 
flected in Figure where gradients are seen in the ratio of 
red to blue clusters as a function of position: the red GCs 
are more centrally concentrated than the blue GCs. Similar ef- 
fects are also seen in other early type galaxies (NGC 1407: 
IFbrbes etal. 20061 NGC 4649 IFbrbes et al. 2()()lat NGC 1399: 
Dirsch et al. 20051 NGC 4636: Dirsch et al. 2005l and others). 
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Since the converse situation is never or rarely seen, this 
must reflect some important difference between blue and red 
clusters. This difference will be related to different formation 
processes, combined with disruption, which will affect popu- 
lations differently depending upon their orbital structure. Less 
radial orbits for the red GCs might explain this, but without 
kinematic data on large samples of GCs nothing conclusive can 
be said about the cause of these differences. 

NGC 2865 shows all signs of a recent merger event: 
a very luminous shell and a KDC with evidence of a re- 
cent, 0.4-1.7 Gyr, st arburst JHauetal. 1999> . and an HI disk 
JSchiminovich 1995i . 

The (V-I) histogram of the GCs of NGC 2865 (FiguregJ is 
more complex than a simple bimodal structure, with a popula- 
tion of very blue, not particularly luminous, GCs near (V-I) = 
0.7. The colour of this population is consistent with an age in 
the range 0.5 - 1 Gyr (Whitmor e et al. 1997> . consistent with 
the nuclear starburst age of 0.4 - 1.7 Gyr f Hau et al. 1999 l. but 
the luminosity of the brightest clusters is much fainter than pre- 
dicted. This could be attributed to the small number of clusters 
in the young population, or else to different physical conditions 
imposing a different Globular Cluster Mass Function. 

An alternative hypothesis is that the structure in the colour- 
magnitude diagrams of the GC systems of NGC 2865 and 
NGC 7626 is entirely due to metallicity variations. In the case 
of NGC 2865 this would require a population of old clusters 
with {V - 1) ~ 0.7, which is too blue for an old population 
at any metallicity (e.g. p:.ee & Carney 200H . In NGC 7626 an 
old population of intermediate metallicity could fill in the dip 
in the CM diagram, but it would have to have a very unusual 
luminosity function to produce the numbers of bright clusters 
that we see. 

7.4. Ages and minor mergers. 

In this Section we compare ages derived for the stellar popula- 
tion of the galaxy; for the shell-forming merger events; and for 
the young GCs that we argue are present in NGC 2865 and 
NGC 7626. Stellar ages derived by comparison with Single 
Stellar Population (SSP) models are uncertain for a number 
of reasons. First, galaxies are clearly not SSPs, and merger 
remnants in particular will have at least three episodes of star 
formation, corresponding to the two progenitors and to the 
merger induced star formation event. Second, the degeneracy 
between age and metallicity, combined with uncertainty in 
isochrone models, and assumptions made about other param- 
eters of the stellar population, such as the mass function, ren- 
der SSP ages very uncertain ( |Poggianti et al. 2001) . Third SSP 
ages are derived from nuclear spectra only, which can be unrep- 
resentative of the galaxy as a whole (e.g. Proctor et al. 2005 1. 
Nevertheless SSP ages have been computed for four of our 
sample by Denicolo et al. (2005), for NGC 1 344 by Kuntschner 
et al (2002), and for NGC 2865 a starburst age has been derived 
by Hau et al. (1999), and these are listed in Table 5. 

Schweizer & Seitzer (1992) determine a Fine Structure 
Index, based upon the galaxy morphology, and derive an em- 
pirical correlation between this and the stellar age. The Fine 



Structure Parameter measures the age since the galaxy as a 
structure was built up, and this is different from the ages of the 
stars in the galaxy. However it determines a combination of the 
age of and the magnitude of a merger event, so there is a degen- 
eracy here as well. Schweizer & Seitzer quote Fine Structure 
ages for three of our sample which are listed in column 3 of 
Table 5. A similar approach was carried out by Michard & 
Prugniel (2004). They divided their peculiar elliptical galaxy 
sample into a normal (NP), reddish, sample, with no signs of 
a young stellar population and a bluish sample (YP) with evi- 
dence for a younger stellar population mixed with an old one. 
They list five of our galaxies and their results are listed in the 
fourth column of Table 5. 

Dynamical ages for phase-wrapped, type I shell systems 
such as NGC 1344, NGC 3923 and NGC 5982 can be derived 
from shell radii and spacings of the outer shells, where dy- 
namical friction and tidal effects are unimportant in determin- 
ing the particle distribution fNulsen 1989 1. We used Nulsen's 
equation 88 to derive the dynamical ages. Mass estimates were 
taken from the literature (NGC 1344: Teodorescu et al. 20051 
NGC 3923: lFuk azawa et al. 2006 ; fundamental plane masses: 
eq. 7 of Ivan D okkum & Stanford 2003 ). The calculated ages 
are shown in Table 5 and are of the order of several 100 Myrs. It 
is clear that these young ages are not reflected in the GC popu- 
lations for these galaxies, implying that the shell forming event 
did not in these cases give rise to a substantial GC population. 

Our sample consists of minor merger remnants, and only 
if it is possible to separate the old and merger related popula- 
tions, as in the case of NGC 2865, can we see a correspondence 
between the different age estimators. More detailed studies of 
the stellar populations of all of the other galaxies would be 
valuable, as would better theoretical estimates of the dynam- 
ical ages of the more complicated Prieur (1990) Type II and III 
shell systems. 

8. Conclusions 

The properties of the globular cluster systems of six shell galax- 
ies were analysed. For NGC 2865 and NGC 7626 we observe 
anomalous features in the (V-I) histograms, in addition to the 
bimodal structure which is normal for ellipticals. The features 
represent excesses at intermediate colour in NGC 7626, and at 
very blue colours in NGC 2865, consistent with a small popu- 
lation of GCs, of age 2-5 Gyr and 0.5-1 Gyr respectively, possi- 
bly formed in the merger event which created the shells. In each 
case the young population is dominated by the much larger, old, 
bimodal distribution. 

The data for two galaxies (NGC 1344 and NGC 3923) al- 
low the determination of their globular cluster luminosity func- 
tion distances. Fitting Gaussians to their GCLF give distances 
moduli of 31.72 and 31.50 for NGC 1344 and NGC 3923 re- 
spectively. 

The properties of NGC 3923 and NGC 5982 are very sim- 
ilar. Their bimodal V-I distributions and radial density profiles 
of blue and red clusters are typical for old GC systems in elUp- 
ticals. NGC 474 and NGC 1344 show one single blue peak and 
very shallow red peaks in their V-I histograms. These properties 
are unusual for bright ellipticals CKundu & Whitmore 200 lal 
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[Peng et al. 2006^ and may indicate less early mergers in their 
formation history. 

NGC 3923 and NGC 7626 have higher specific frequencies 
(respectively 5.6 and 3.9 within 50kpc) than normal for galax- 
ies located in a low density environment. The luminosity of the 
X-ray halo detected in the former galaxy is probably not suffi- 
cient to explain its high Sn proposed by McLaughUn (1999)| 
The Sn of the other galaxies have values of around 2 and are 
typical for galaxies located in a low density environment. 

Although for some of these six galaxies new GCs may have 
been formed recently, the general properties (like V-1 distribu- 
tions and flattening of GC density profile) of the globular clus- 
ter systems of these shell galaxies do not deviate systematically 
from 'normal' elliptical galaxies. 

In a future paper (Sikkema et al., in preparation) we will in- 
vestigate the morphology and stellar populations of the diffuse 
galaxy, in particular in the shells. 
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